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Abstract. - The Modified Fluctuation Dissipation Theorem (MFDT) proposed by G. Verley et al. 
(EPL 93, 10002, (2011)) for non equilibrium transient states is experimentally studied. We apply 
MFDT to the transient relaxation dynamics of the director of a liquid crystal after a quench close 
to the critical point of the Freedericksz transition (Ftr), which has several properties of a second 
order phase transition driven by an electric field. Although the standard Fluctuation Dissipation 
Theorem (FDT) is not satisfied, because the system is strongly out of equilibrium, the MFDT is 
perfectly verified during the transient in a system which is only partially described by Landau- 
Ginzburg (LG) equation, to which our observation are compared. The results can be useful in the 
study of material aging. 


C3‘ 


T3 

G 

O 

O 


> 

'xl- 

'yI- 

y|- 

m 

o 

o 

i n 


X 


After a sudden change of a thermodynamic parameter, 
such as temperature, volume and pressure, several systems 
and materials may present an extremely slow relaxation 
towards equilibrium. During this slow relaxation, usually 
tailed aging, these systems remain out-of equilibrium for 
a very long time, their properties are slowly evolving and 
equilibrium relations are not necessarily satisfied during 
aging. Typical and widely studied examples of this phe¬ 
nomenon are glasses and colloids where many questions 
on their relaxation dynamics still remain open [1,2]. Thus 
in order to understand the minimal ingredients for aging, 
slow relaxations have been studied theoretically in sec¬ 
ond order phase transitions when the system is rapidly 
quenched from an initial value of the control parameter 
to the critical point [3-6]. Because of the critical slow¬ 
ing down and the divergency of the correlation length the 
relaxation dynamics of the critical model shares several 
features of the aging of more complex materials. One of 
the questions analyzed in this models is the validity of 
the Fluctuation Dissipation Theorem (FDT) during the 
out of equilibrium relaxation [8-11]. In equilibrium, FDT 
imposes a relationship between the response of the system 
to a small external perturbation and the correlation of the 
spontaneous thermal fluctuations. When the system is out 
of equilibrium FDT does not necessarily hold and it has 


been generalized as 

k B T X(t,t w ) x(t,t w ) 


C{t, t) - C(t, t w ) (1) 


where k B is the Boltzmann constant, T the bath temper¬ 
ature, C(t,t w ) =< 0(t)0(t w ) > - < 0(t) >< 0(t w ) > 
(< . > stands for ensemble-average) the correlation func¬ 
tion of the observable 0(t). The function x(i, t w ) is the 
response to a small step perturbation, of the conjugated 
variable h of O, applied at time t w < t : 




< 0(t ) > h - < 0(t) > 0 
h 


I h —>-0 


( 2 ) 


where < 0(1) >h and < 0(t) > a denote respectively the 
mean perturbed and unperturbed time evolution. The 
function X(t,t w ) is equal 1 in equilibrium whereas in out- 
equilibrium it measures the amount of the FDT violation 
and it has been used in some cases to define an effective 
temperature T e ff(t,t w ) = X(t,t w )T. 

The above mentioned models of the quench at critical 
points allows a precise analysis of the pertinence of this 
definition of T e ff [8]. In spite of the theoretical interest of 
these models only one experiment has been performed on 
the slow relaxation dynamics after a quench at the critical 
point [7]. The role of this letter is to experimentally ana¬ 
lyze the theoretically predictions in a real system affected 
by finite size effects and unavoidable imperfections. We 
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also analyze another important aspects of the FDT in out 
of equilibrium systems. Indeed several generalizations of 
FDT has been proposed [9-17] but almost all of them can 
be applied to non-equilibrium steady states [18,19] and 
are not useful for the transient time evolution which fol¬ 
lows the quench at the critical point. As far as we know 
there are only two formulations of FDT [20,21], which are 
useful for these transient states and the second purpose 
of this letter is to discuss the application of the Modified 
FDT (MFDT) of ref. [20]. 

Before describing the experimental set-up we summa¬ 
rize briefly the formulation of the MFDT of ref. [20] . Let 
us consider the relaxation dynamics of a system, which 
has been submitted at time t = 0 to a sudden change of 
its control parameters. At time t after the quench, this 
relaxation is characterized by the variable x(t), by the ob¬ 
servable 0(x(t)) and by the probability density function 
n(x(t), h(t)). Here h(t) is an external control parameter 
which is used to perturb the dynamics. We define a pseu¬ 
dopotential 'F(t, h) = — ln[7r(a;(t), h)] and an observable 
B{t) = —dh4>{t, h) Ih-un with h^O constant for t > 0 and 
h = 0 for t < 0. The MFDT reads: 

= < B(t)0(t) > - < B(t w )0{t) > (3) 

Eq.3 defines the response function of 0(t) to a 

step perturbation of h applied at t w , with 0 < t w < t. 
Notice that in this case h can be any parameter of the 
system and it does not need to be the conjugated variable 
of 0(t). In this letter we will analyze how MFDT (eq.3) 
can be applied to the experimental data using the quench 
at the critical point in a non ideal system. 



Fig. 1: Phase diagram of the Freedericksz Transition in 5CB. 
a) Dependence on e of 6 2 =< 6 > 2 (blue dots) and of the 
variance ( dashed red line) of 6 multiplied by 4 10® to be on 
the same scale. The solution of the LG equation (p = 0) is 
the straight black line, b) t 0 /t is plotted as a function of 
t. The prediction of LG (p = 0) is the black straight line. 
The measured relaxation time deviates significantly from the 
predicted one even for values of e where the stationary solution 
of LG, plotted in a) seems to reproduce the data. 

The experimental system where we study these proper¬ 
ties and the MFDT is the Freedericksz Transition (FrTr) 
in a nematic liquid crystal (LC) submitted to an external 
electric field E. Specificlly in our experiment we use the 


5CB (p-pentyl-cyanobiphenyl, 5CB, produced by Merck). 
The experimental apparatus has been already described 
[22, 23] and we summarize here only the main features. 
The LC is confined between two glass plates, separated by 
a distance L =13.5 pm. The surfaces in contact with LC 
molecules are coated by ITO to apply an electrical field. 
Then, a polymer layer (rubbed PVA) is deposited to in¬ 
sure a strong anchoring of the 5CB molecules in a direction 
parallel to the plates. In the absence of any external field, 
the molecules in the cell align parallel to those anchored at 
the surfaces. Applying a voltage difference U between the 
electrodes, the liquid crystal is submitted to an electrical 
field perpendicular to the plates. To avoid polarization, 
the applied voltage is modulated at a frequency f = 10 
kHz [U = v / 2L f 0 cos(27r/t)]. When U 0 exceeds a thresh¬ 
old value U c , the planar states becomes unstable and the 
molecules rotate to align with the electrical field. To quan¬ 
tify the transition we measure the spatially averaged align¬ 
ment of the molecules determined by the angle 6 between 
the molecule director and the surface. Such measurement 
relies upon the anisotropic properties of the nematic. This 
optical anisotropy can be precisely measured using a very 
sensitive polarization interferometer [22] which gives a sig¬ 
nal <p tx 9 2 . At U ~ U c the dynamics is usually described 
by a Landau-Ginzburg (LG) equation although as pointed 
out in ref. [23] this is a very crude approximation, which 
has several drawbacks. In a very first approximation the 
dynamics of the mean relaxation 0o(t) = (0(f)) is ruled by 
the following Ginzburg-Landau equation : 

t o 0q = ed 0 - | ( 0 q - p 3 ) ( 4 ) 


where e = (U 2 — U 2 )/U 2 is the reduced control parameter, 
the t 0 is a characteristic time of the LC and a a parameter 
which depends on the elastic and electric anisotropy of 
the LC. For 5CB a = 3.36 and t 0 = 2.4s for the cell 
thickness L = 13.5 pm. The residual angle p ~ 0.1 at 
e = 0 comes from cell assembling and preparation and has 
been discussed in ref. [23]. Furthermore r Q is not strictly 
constant but it slightly depends on 9 2 . The dimensional 
equation for the fluctuations 86{t) = 6(t) — 9 0 (t) is 


r )AL59 
with K 


K[(e-'^a9l) 89 + Se9 0 \+ v (5) 

TT 2 klA 

— (6) 


where A is the laser cross section, k\ is one of the elas¬ 
tic constant of LC and rj is a delta correlated thermal 
noise such that < r/(t)r](t') >= kBT(jAL)5(t — t'). The 
term with Se takes into account that during the measure 
of the response we have to perturb the value of e by ap¬ 
plying a short pulse of duration r p and amplitude Se. The 
two eqs.4,5 describe in principle the dynamics of the mean 
deflection 9 a and of the fluctuations 56. However there 
are several discrepancies with the experimental data which 
are widely discussed in ref. [23]. We summarize here the 
most important, which are useful for the discussion. The 
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phase diagram of FrTr in 5CB and the relaxation times 
are plotted in figs. 1 a),b) respectively. The solution of 
LG 91 ~ 2e/a (/i = 0 in eq.4) reproduces the station¬ 
ary experimental data for e < 0.6. Instead the measured 
relaxation time (fig. 1 b) deviates significantly from the 
predicted one even for values of e where the stationary 
solution of LG, (fig. 1 a) seems to reproduce the data. 
The characteristic time increases but it does not diverge 
because of /r ^ 0 (see fig.l and ref. [23]). In fig. 1 a) the 
variance Ug of 6 is plotted too. From eq.5 this variance is 
= fcsT/(A'(e— | a9g)) which does not diverge at e = 0 
because fi ^ 0. 

Thus although eq. 4 and eq.5 are only a rough approx¬ 
imation of the FrTr dynamics, especially at e > 0.1, we 
use them to fix the framework, and because they are very 
close to the theoretical mean field approach to the quench 
at critical point discussed in ref. [6] . Thus it is interesting 
to check the analogies and differences with respect to the 
general theory. 

The quench is performed by commuting e from an initial 
value Ci to an ef ~ 0 at t = 0. As an example we show 
in fig. 2 the time evolution of 0 o for a quench from et = 
0.25 to e/ = 0.01. The system is relaxing from its initial 
equilibrium value towards the new one. We describe here 
the time evolution of the statistical properties and we will 
discuss at the end the dependence on the initial and final 
e values. The mean values of the statistical properties are 
obtained by repeating the quench at least 3000 times. 



Fig. 2: Quench close to the critical point from a = 0.25 to 
e/ = 0.01. a) Time evolution of the order parameter 9 2 before 
and after the quench preformed at t = 0. b) Time evolution of 
the variance as a function of time. The variance and the mean 
9 2 (t) has been obtained by performing 3000 quenches and then 
making an ensemble average on the quenches at each time. 

The time evolutions of 9^(t) and of the variance <Jg(t) 
are shown in figs. 2a) and b). We see that both quantities 
relax from the initial to the final equilibrium values ,which 
are 9\ ~ 2e/a and crg{t) = ksT/ (A'(e— § a0^)), where 9 2 is 
the equilibrium value, which is not exactly 2 e/a because 
of the presence of the imperfect bifurcation fi / 0 (see 
fig.l and ref. [23]). We see that the fluctuation amplitude 
increases when approaching the critical point. 

In Fig. 3 we plot C(t, t w ) =< S9(t)S9(t w ) > as a function 
of t — t w at various t with t > tw > 0. We see that C(t, t w ) 
develops very long decays when t is increased. In order to 



Fig. 3: The correlations functions (a) and the integrated re¬ 
sponses (b) (computed at various fixed times t and 0 <t w <t 
during the relaxation after the quench) are plotted as a func¬ 
tion of t — tw 


study FDT we need to measure the response by perturb¬ 
ing the systems with a pulse of amplitude 6e = 0.4 and 
Tp = 1 ms at time t w . As an example, in fig.4a) we plot the 
time evolution perturbed at t w = 5 s and in fig.4b) the time 
evolution of the difference < A 9{t) >=< it) > —9 0 (t) 
between the perturbed 9sit) and the unperturbed 9 0 {t). 
As it can be seen in eq.5 the amplitude of the pertur¬ 
bation is K5e{t w )9{t w )T p . Thus the impulse response 
function is R{t,t w ) =< A 9(t) > /(< 5e(t w )9(t w )T p ) for 
t w < t. We repeat the experiments N p times by send¬ 
ing at each quench a pulse at a different time t Wt i with 
[^io,i = 0,., t w ,Np = 20s]. Then the integrated re¬ 

sponse is 

iVt-l 

X'(L t w ^ m ) — ^ t w ,i-i-l) (tyj,(i-i~l) t'UJ,i) (7) 

i—m 

such that t = t w ,N t and \(t, t) = 0. 

The measured x(f, t w ) is plotted as a function of t — t w for 
various t in fig.3.b. 

To check the validity of the standard FDT, we plot, in 
fig. 5, x{t>tw)kBT/C(t,t) as a function of C(t,t w )/C(t,t) 
at various fixed t with t w varying in the interval 0 < t w < 
t. In this plot FDT is a straight line of slope -1. We see 
that for t relatively short, compared to r, the FDT is not 
satisfied. In fig. 5 we also plot the prediction of ref. [6] 
for a quench done at e = 0 in a Landau-Ginzburg (LG) 
equation. We see that for short time the behavior is quite 
different from that of the LG equation confirming that the 
dynamics is not very well described by this equation. The 
behavior at long time is instead related to the fact that 
the quench is not performed exactly at e/ =0. 

We now apply the MFDT to these data. In order to do 
that one has to considers that 69 has a Gaussian distri¬ 
bution whose variance is plotted in fig. 2. As observable 
in eq.3 we use 0(x(t)) = x{t) = 69. Following the for¬ 
mulation of the MFDT one has to consider the dynamics 
of 'F(t) when a small perturbation h is applied at t = 0, 
therefore by the definitions of x(M™) ( ec l-7) an d of 0(t), 
we get < 69{t) >h= x(^> 0) h because < 69{t) > 0 = 0. 
Thus at h ^ 0 (switched on at t = 0) the probability 
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Fig. 4: a) The time evolution has been perturbed at time t w 
by a short pulse of amplitude St = 0.4 and duration 1ms. b) 
The response A 9 to the delta perturbation. 



Fig. 5: FDT plot. The function x(t, t w )kBT/C(t, t) is plotted 
as a function of C(t,t w )/C(t,t) at various fixed times t and 
0 < tw < t after the quench. In equilibrium this plot is a 
the straight line of slope —1. We see that at short times t 
the curves strongly deviate from the equilibrium position. The 
equilibrium FDT is recovered only for very large t. The red 
dashed straight line is the prediction [6] for a quench done at 
£ = 0 in a Landau-Ginzburg equation. 


density function for of 89 around the mean is 

n(S9(t),h) = V / l/(27rg 2 (f))exp(^ ^ 61 ^ 2 ^ ^ )> ( 8 ) 

where we assume that if h is small enough then the depen¬ 
dence of crg(t) on h can be neglected. 1 . Therefore from 
the expression of 7r (8d(t),h), the definition of ^(89(t), h) 
and of B(t) one finds: B(t) = 89(t)x{t, 0)/ crg(t ). Thus eq. 
3 for this particular choice of variables becomes : 

x(b _ C(f) tw'jxitw i 0) 

All the quantities in eq.9 have been already measured. 
Thus in fig. 6 we plot x(t>t w )/x(t,0), the left term of eq.9, 
as a function of C(t, t w )x(t w , 0)/(x(f, 0 )<Jg(t w )) for various 
fixed t and 0 < t w < t. We see that all the data points are 
aligned on a straight line of slope —1 as predicted by eq.9. 

1 This has been verified experimentally and for the dynamic de¬ 
scribed by eq.6 



Fig. 6: In order to verify MFDT the left hand side 

of eq.9, i.e. x(f, t w )/x{t, 0), is plotted as a function of 
C(t,t w )x{tw,0)/x(t,0)&e(t w ). All the data collapse on the 
straight line of slope —1 showing that the MFDT prediction 
for transient is perfectly verified for any t. Notice that the re¬ 
sponse and the correlations are the same than those used in 
fig. 5, but now they have been normalized as prescribed by the 
MFDT, i.e. eq.9. 



Fig. 7: MFDT recomputed using in eq.9 the x!{tR w ) defined 
in the text. The left hand side of eq..9 is plotted as a function 
of the right hand side. All the data collapse on the straight line 
of slope —1 showing that the MFDT prediction for transient is 
perfectly verified for any t in this case too. 


We clearly see that, in contrast to fig. 5 where the stan¬ 
dard formulation of FDT is recovered only for very large 
t , MFDT is verified for all times. As pointed out there is 
no need in MFDT to use for h the conjugated variable of 
0{t). We can use simply h = <5e. In such a case we define 
the response function as R!{t,t w ) =< A 9 > /{8e t p ) and 
the x'(titw) is obtained by inserting R'{t,t w ) in eq.7. The 
MFDT computed using in eq.9 x'(t, t w ) instead of x(b t w ) 
is checked in fig. 7 where the left hand side of eq.9 is plot¬ 
ted as a function of the right hand side. We see that the 
MFDT is verified in this case too . 

All the data presented in this paper correspond to a 
quench from ti ~ 0.25 to £/ ~ 0.0, however the main 
statistical features, here described, are independent on the 
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starting and final points. The final point influences the 
duration of the out of equilibrium state, which depends on 
the distance from the critical point. The small difference 
with the results in ref. [7] is due to a slightly non linear 
response in that reference. 

As a conclusion in this letter we have applied to a 
quench a the critical point of Freedericksz transition (Ftr), 
the Modified Fluctuation Dissipation Theorem for tran¬ 
sient, proposed in ref [20]. We find that although the 
equilibrium FDT is strongly violated the GFDT is very 
well satisfied, independently of the chosen response. It 
is interesting to point out that the result is interesting 
because although the system presents several differences 
with respect to the LG equation, it is affected by finite 
size effects and the quench is not performed exactly at the 
critical point the dynamics still presents features at the 
critical quenching. 

We acknowledge useful discussion with G. Verley and 
D. Lacoste. This work has been supported by the ERC 
contract OUTEFLUCOP. 
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